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ABSTRACT 

We present the Pan-STARRSI discovery of PSl-lOafx, a unique hydrogen-deficient superluminous 
supernova (SLSN) at redshift z = 1.388. The hght curve peaked at zpi — 21.7 mag, making PSl-lOafx 
comparable to the most luminous known SNe, with Mu — —22.3 mag. Our extensive optical and near- 
infrared observations indicate that the bolometric light curve of PSl-lOafx rose on the unusually fast 
timescale of ^12 d to the extraordinary peak luminosity of 4.1x10''^ erg s~^ (Mboi = —22.8 mag) and 
subsequently faded rapidly. Equally important, the SED is unusually red for a SLSN, with a color 
temperature of ^^6800 K near maximum light, in contrast to previous hydrogen-poor SLSNe, which 
are bright in the ultraviolet (UV) . The spectra more closely resemble those of a normal SN Ic than any 
known SLSN, with a photospheric velocity of '^ll, 000 km s~^ and evidence for line blanketing in the 
rest-frame UV. Despite the fast rise, these parameters imply a very large emitting radius (>5xl0^^ 
cm). We demonstrate that no existing theoretical model can satisfactorily explain this combination 
of properties: (i) a nickel-powered light curve cannot match the combination of high peak luminosity 
with the fast timescale; (ii) models powered by the spindown energy of a rapidly-rotating magnetar 
predict significantly hotter and faster ejecta; and (iii) models invoking shock breakout through a dense 
circumstellar medium cannot explain the observed spectra or color evolution. The host galaxy is well 
detected in pre-explosion imaging with a luminosity near L*, a star formation rate of ~15 MQyr~^, 
and is fairly massive (~2xl0^° M0), with a stellar population age of ~10^ yr, also in contrast to the 
young dwarf hosts of known hydrogen-poor SLSNe. PSl-lOafx is distinct from known examples of 
SLSNe in its spectra, colors, light-curve shape, and host galaxy properties, suggesting that it resulted 
from a different channel than other hydrogen-poor SLSNe. 

Subject headings: supernovae: individual (PSl-lOafx) — surveys (Pan-STARRSI) 



1. INTRODUCTION 

A small fraction of massive stars end their lives with 
spectacular explosions one or two orders of magnitude 
more luminous than normal supernovae (SNe). After the 
initial puzzling discoveries of the luminous SNe 2005ap 
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(IQuimbv et all l2007l) and 2006gy (|Ofek et al.l [20071: 
ISmith et all 120071) . modern wide-field surveys over the 
past decade began to uncover these superluminous SNe 
(SLSNe) in greater numbers. The energy scales involved 
in these explosions challenge our understanding of con- 
ventional SN explosions. Normal SNe resulting from iron 
core collapse have characteristic energy scales of ~10^^ 
erg of kinetic energy and ^10'*^ erg emitted as optical 
radiation (^lO'*^ erg s~^ for ^10'' s). The SLSNe are far 
off this scale — they peak at optical luminosities of up 
to ~4xl0^^ erg s~ ^ (Quimby ct al. 2011; Chomiu k et al.l 
[20Tlt iMiUer et all [20091: iGezari et all 120091) and em it a 
total of up to 4x10^1 erg optically ([Rest et al.ll201lD . 

This large energy scale motivates the question of what 
physics powers these SNe, and how to accommodate 
these objects within the standard understanding of mas- 
sive star evolution. Theorists have proposed a number of 
exotic powe r sources, including t he pair instability mech- 
anism (e.g., IBarkat et ahl 119671: IRakavv fc Shavivlll967[) 
and reprocessed spi ndown energy r el eased by a newly 
forme d magnetar ([Wooslevl [20101 : [Kasen fc BildsteiJ 
[2010| ). Another possibility is interactio n with a 
dense circumstellar medium (CSM ) (Chevalier fc Irwin 



2011 1 : Moriya & T ominaeal [2011 IGinzburg fc Balbere 
20121 : FChatzopoulos et al.ll2012[ ). requiring extreme CSM 



mas ses and densities wh ose origin remains unexplained 
(see IWooslev et all[2007[ for one possibility) . All of these 
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models require additional ingredients beyond the normal 
st ellar evolutiona ry processes. 

iGal-YarnI ()2012[ ) has attempted to impose order on the 
menagerie of objects achieving sufficient peak luminosi- 
ties to be classified as SLSNe (A/ < — 2f mag is a typical 
requirement) by sorting them into three categories. All 
of the hydrogen-rich objects were classified as SLSNe-II 
and all exhibit signs of being powered by dense GSM 
interaction, with t he possible except ion of SN 2008es 
(jMiller et al.ll2009HGezari et al.ll2009[ ). He split the ob- 
jects lacking hydrogen into two classes, the rare SLSNe- 
R that have slow photometric decline rates consistent 
with being powered by the radioactive decay of a very 
large synthesized mass of ^^Ni, and the relatively ho- 
mogeneous class of SLSNe-I, whose power source is still 
mysterious. A few caveats have been raised. The SLSNe- 
R are interpreted to be the results of pair-instability 
SNe. However, existing models for the pair instabil- 
ity process prefer extremely low metallicity, and may be 
in conflict with the observed spectrum an d spectral en- 
ergy d istribution (SED) of SLSNe-R fe.g.. iDessart et al.) 
I2012b( ). Also, it is not clear how homogeneous the 
SLSNe-I class really is. Although the spectra of most ap- 
pear to be similar to those of SN 2 005ap and SCP06F6 
(jQuimbv et "alll2007l: iBarbarv et al.|[2009i ). the rise times 
and peak lu minosities of p ublish ed objects vary by fac- 
tors of ^ 5 (iQuimby et al.l 120 07!: iPastorello et all 120101: 
Ghomiuk et al. 2011; Lclouda s et al.l 
All SLSNc-I to date have 
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had hot spectra and been bright in the rest-frame near- 
ultraviolet (NUV) relative to normal SN SEDs. 

In this paper, we present the discovery of PSl-lOafx, 
an extreme SLSN at redshift z = 1.388 that does not 
fit into this classification scheme and is distinct from all 
previous SLSNe. The peak luminosity is comparable to 
the highest known and the rise time is the fastest mea- 
sured. The spectra show no evidence for hydrogen and 
lack any analog in the existing sample of SLSNe. Instead, 
they most closely resemble those of line-blanketed nor- 
mal SNe Ic. In Section 2, we present the suite of optical 
and near-infrared (NIR) observations. The host galaxy 
is described in Section 3. We compare our observations 
of PSl-lOafx to known SNe in Section 4. In Section 5, we 
construct the SED and bolometric light curve. We then 
compare PSl-lOafx to existing SLSN models in Section 
6. All calculations in this paper assume a flat ACDM 
cosmolog y with Hn—7A km s~ ^ Mpc~^, r2m=0.2 7, and 
r2A=0.73 (|Komatsu et al.l[20TTI : iRiess et al.ll20Tl[ ). 

2. OBSERVATIONS 
2.1. Pan-STARRSl 

The Pan-STARRSl (PSl) telescope has a 1.8 m di- 
ameter primary mirror that images a fleld with a di- 
ameter of 3.3° ()Hodapp et all I2004D onto a total of 
sixty 4 800 x 4800 p ixel detectors, with a pixel scale of 
0.258" (|Tonrv fc On aka 2009). A more complete descrip- 
tion of t he PSl sys tem, hardware and software, is pro- 
vided bv lKaiser et'a l. (2010). 

The PSl observations are obtained through a set of 
five broadband filters, designated as gpi, rpi, ipi, 2;pi, 
and 2/pi. Although the filter system for PSl has much 
in common with that used in previous surveys, such as 
the Sloan Digital Sky Survey (SDSS; ,Ahn et al.,,2012, ). 



there are differences. Most important for this work, the 
zpi filter is cut off at 9300 A, giving it a different re- 
sponse than the detector response defined zsdss^ and 
SDSS has no corresponding ?/pi filter. Fur ther informa- 
tion o n the passband shapes is described bv lStubbs et al.l 
(|2010f) . Photometry is in the "natural" PSl system, 
m = — 2.51og(/y) + to', with a single zeropoint adjust- 
ment to' mad e in each band to c onform to the AB mag- 
nitude scale (jTonrv et al.l 120121 ). Photometry from all 
other sources presented in this paper is also on the AB 
scale. PSl magnitudes are interpreted as being at the 
top of the atmosphere, with 1.2 airmasses of atmospheric 
attenuation being included in the system response func- 
tion. 

The PSl Medium Deep Survey (MDS) consists of 10 
fields across the sky that are observed nightly when in 
season (■^S months per year) with a typical cadence of 
3 d between observations in gpirpiipizpi in dark and 
gray time, while ypi is used near full moon. PSl data 
are processed through the Image Processing Pipeline 
(IPP; Magnier 2006) on a computer cluster at the Maui 
High Performance Computer Center. The pipeline runs 
the images through a succession of stages, including 
fiat-fielding ("de-trending"), a flux-conserving warping 
to a sky-based image plane, masking and artifact re- 
moval, and object detection and photometry. Tran- 
sient detection using IPP photometry is carried out at 
Queens University Belfast. Independently, difference im- 
ages are produced from the sta cked nightly MD S im- 
ages by the photpipe pipeline (iRest et al.(l2005D run- 
ning on the Odyssey computer cluster at Harvard Uni- 
versity. The discovery and data presented here are from 
the photpipe analysis. 

PSl-lOafx was first detected in MDS imaging on 2010 
August 31.35 (UT dates are used throughout this paper) 
at a position of a = 22hll'°24!162, 5 = +00°09'43'.'49 
(J2000), with an uncertainty of O'.'l in each coordinate. 
The strong detections in ipi and zpi combined with non- 
detections in gpi and rpi over the next few nights im- 
mediately garnered our attention. The unusual colors 
are evident in the color images of the field presented in 
Figure [TJ 

We constructed deep template stacks from all pre- 
explosion images and subtracted them from the PSl ob- 
servations using photpipe (Rest et al. 2005). Details of 
the photometry and generation of PSl SN light curves 
will be given by Rest et al. (2013, in prep.) and Scol- 
nic et al. (2013, in prep.). The typical spacing be- 
tween MDS observations in the same filter corresponds 
to only 1.3 d in the rest frame of PSl-lOafx, so in some 
cases we have co-added the photometry from adjacent 
observations to increase the significance of marginal de- 
tections or the depth of the non-detections. The final 
PSl-lOafx photometry, after corr ection for E(B — V) = 
0.05 mag of Galactic ext inction (jSchlaflv fc Finkbeinei] 
l201lHSchlegel et al.l[T99l . is given in Tableland shown 
in Figure [2] 

We fit polynomials to the zpi data points within 15 
rest-frame days of peak to determine the time of maxi- 
mum light, for which we adopt a Modified Julian Date 
(MJD) of 55457.0 (=2010 September 18.0). Ah phases 
referred to subsequently are in rest-frame days referenced 
from this date. 
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Figure 1. Three-color gpirpiipi images of the field of PSl-lOafx, showing {left) a deep stack of prc-explosion imaging witli the host 
galaxy marked; (center) images taken near maximum light; and {right) difference images of the templates subtracted from the observations. 
The color scales are similar in each panel. The unusually red color of PSl-lOafx compared to other faint objects in the field is apparent, 
as it is only strongly detected in ipi, with non-detections in gpi and rpi. 
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Figure 2. Multicolor optical light curves of PSI-lOafx. Most of 
the data were obtained with PSl, but points marked r', i' , or z' 
come from other telescopes (see Table [TJ. The r'i'z' photometry 
points near days —3 and -|-3 were taken contemporaneously with 
the optical spectra shown in Figure |4] Upper limits marked with 
arrows correspond to 3cr. 



2.2. Other Photometry 

In addition to the PSl observations, we obtained two 
epochs of multicolor photometry using the Gemini Multi- 
Object Spectrographs (CMOS; Hook et al. 2004) on the 
8-m Gemini- North and South telescopes and one epoch of 
imaging using the Ina mori-Magellan Areal Camera and 
Spectrograph flMACS: [D"ressler et al.ll2006D on the 6.5- 
m Magellan Baade telescope. The images were processed 
using standard tasks and then archival fringe frames were 
subtracted from the GMOS images using the gemini 
IRAFI3 package. We used several SDSS stars in the 
field to calibrate the Gemini images, while the IMACS 
zeropoints were checked with observations of standard 
star fields obtained the same night. We subtracted the 
deep PSl templates in the corresponding filter from each 
imag e using the ISIS software package (Alard & Lupto^ 
119981 ) to correct for host contamination. The final pho- 
tometry is listed in Table [1] The r' and i' magni- 
tudes agree well with the PSl observations at similar 
epochs. However, the z' observations exhibit an offset 
from the zpi light curve due to the filter response dif- 
ferences noted above (GMOS and IMACS are closer to 
SDSS). We therefore add 0.2 mag to the z' photometry 
points for consistency with zpi whenever we refer to the 
combined z-band light curve. Finally, we obtained late- 
time observations of the host galaxy in g' and r' on 2011 
October 21.1 (phase -f 167 d) using the Low Dispersion 
Survey Spectrograph-3 (LDSS3) on the 6.5-m Magellan 
Clay telescope (Table [2]). 

The high redshift and red optical colors of PSl-lOafx 
indicate that NIR photometry is required to constrain 
the SED and bolometric luminosity. We received direc- 
tor's discretionary (DD) time at Gemini- North to obtain 
YJHK photometry near maximum light using the Near 
Infra Red Imager and Spectrometer (NIRI: [Hodapp et"ari 
|2003[ ). The observations were acquired on 2010 Septem- 
ber 24.4, at a phase of -1-2.7 d. Additional observations 
were obtained over the subsequent month in JHKg with 
the MMT & Magellan Infrared Spectrograph (MMIRS; 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1 

PSl-lOafx Photometry 



MJD 


Epoch^ 


Filter 


Magnitude'' 


Error 


Instrument 




(d) 




(AB) 






55440.5 


—6.9 


9P1 


>23.50 




PSl 


55445.0 


—5.0 


9P1 


>23.91 




PSl 


CC/ICn QC 




9P1 


^Zo. / / 






cr c /I c c /I c 

oo4o5.4 


—0.7 


9P1 


>23.22 




PSl 


5547U.O 


5.6 


9pi 


>23.68 




PSl 


r r /I Qn '7C 


9.9 


9pi 


>23.85 




Pol 


55439.0'- 


—7.5 


rpi 


>23.13 




rrrn 

PSl 


55445.0'^ 


—5.0 


rpi 


>23.47 




PSl 


55450.2 


—2.8 


r' 


24.18 


0.30 


GMOb-N 


55450.9 


—2.5 


^•pi 


23.87 


0.41 


PSl 


55455. 4"^ 


—0.7 


rpi 


>23.42 




PSl 


55465.0 


3.4 


r' 


23.73 


0.20 


GMOS-S 


55470.3'^ 


5.6 


rpi 


>23.63 




PSl 


55480.8 


9.9 


^•pi 


>23.55 




PSl 


55425.0*- 


— 13.4 


«P1 


>24.21 




PSl 


55438.5"^ 


—7.7 


jpi 


>23.84 




PSl 


55441.5 


-6.5 


ipi 


23.10 


0.16 


PSl 


55444.4 


-5.3 


jpi 


22.57 


0.12 


PSl 


55447.4 


—4.0 


jpi 


22.32 


0.06 


PSl 


55450.4 


—2.8 


ipi 


22.15 


0.11 


PSl 


55450.2 


-2.8 


i' 


22.18 


0.11 


GMOS-N 


0040O.4 


— i.O 


ipi 




U.Uo 


r ol 


55456.4 


—0.3 


jpi 


21.99 


0.05 


PSl 


55465.0 


3.4 


i' 


22.25 


0.16 


GMOS-S 


55465.5 


3.5 


ipi 


22.28 


0.20 


PSl 


55468.3 


4.7 


ipi 


22.39 


0.12 


PSl 


55477.2 


8.5 


ipi 


23.20 


0.19 


PSl 


55480.2 


9.7 


ipi 


23.28 


0.42 


PSl 


55483.2 


11.0 


ipi 


23.45 


0.22 


PSl 


55486.3 


12.2 


ipi 


>22.90 




PSl 


55499.8 


17.9 


ipi 


>23.87 




PSl 


55508.8 


21.7 


ipi 


>24.08 




PSl 


55513.1 


23.5 


■/ 

t 


>25.01 




TTV T A /~l O 

IMAGS 


55426.0*- 


— 13.0 


Zpi 


>23.47 




PSl 


55436.4 


-8.6 


Zpi 


>22.79 




PSl 


55439.3 


-7.4 


Zpi 


23.06 


0.37 


PSl 


55442.5 


-6.1 


zpi 


22.83 


0.37 


PSl 


55445.4 


-4.8 


Zpi 


22.24 


0.15 


PSl 


00440.4 


— O.D 




91 


n 1 9 


p<^ 1 


55450.2 


—2.8 


1 

z 


21.59 


0.13 


GMOb-N 


55454.4 


— 1.1 


Zp-i 


21.69 


0.10 


PSl 


55457.3 


0.1 


Zpi 


21.76 


0.11 


PSl 


004D0.U 


0.4 


z 


91 AA 


n 1 Q 




55466.4 


3.9 


zp\ 


21.82 


0.10 


PSl 


55475.2 


7.6 


zp\ 


22.11 


0.17 


PSl 


55478.2 


8.9 


Zpx 


22.13 


0.17 


PSl 


55481.2 


10.1 


Zpi 


21.98 


0.19 


PSl 


55484.2 


11.4 


zp\ 


22.77 


0.30 


PSl 


55499.2*^ 


17.7 


Zpi 


>23.04 




PSl 


55513.1 


23.5 


z' 


23.83 


0.30 


TTV /f A /^r^ 

IMAGS 


D04DO.4 


Z. t 


ypi 


91 no 


0.10 


DQl 


55492.8*^ 


15.0 


ypi 


>21.95 




PSl 


55523.3 


27.8 


ypi 


>21.82 




PSl 


55463.4 


2.7 


Y 


20.96 


0.10 


NIRI 


55463.3 


2.6 


7 

J 


20.99 


0.17 


MMlHo 


55463.4 


2.7 


J 


21.19 


0.10 


NIRI 


55481.1 


10.1 


J 


22.22 


0.26 


IVliVlllto 


55485.0 


11.7 


J 


22.28 


0.30 


MMIRS 


55496.0 


16.3 


J 


22.19 


0.10 


HAWK-I 


55515.1 


24.3 


J 


22.84 


0.14 


HAWK-I 


55552.0 


39.8 


J 


>21.90 




HAWK-I 


55561.0 


43.6 


J 


>21.90 




HAWK-I 


55463.4 


2.7 


H 


21.20 


0.19 


NIRI 


55485.0 


11.7 


H 


21.05 


0.44 


MMIRS 


55490.1 


13.9 


H 


21.94 


0.24 


HAWK-I 


55515.1 


24.3 


H 


22.28 


0.39 


HAWK-I 


55463.4 


2.7 


K 


21.46 


0.26 


NIRI 


55485.0 


11.7 


Ks 


>21.35 




MMIRS 



^ In rest-frame days relative to maximum light on MJD 55457.0. 
^ Corrected for Galactic reddening. Upper limits are 3ct. 
^ Data point is from a stack of observations taken on multiple nights. 
MJD reported is mean of individual observations in stack. 
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Figures. NIR photometry of PSl-lOafx. The dashed line shows a 
polynomial fit to the well-sampled zpi light curve to guide the eye, 
shifted by 0.7 mag to go through the J photometry near maximum 
light. 



iMcLeod et"aLll2012D on the Magellan Clay Telescope. To 
follow the light curve in the rest-frame optical to late 
times, we obtained a series of JH observations using 
DD time on the 8-m Very Large Telescope (Yepun) with 
the High-Acuity ^yide-fi eld K-band Imager (HAWK-I; 
IKissler-Patig et al.ll2008l). Finally, we used the FourStar 
Infrared Camera (jPersson et al.l I2008D on the Magellan 
Baade telescope to obtain late-time JHKs observations 
of the host galaxy of PSl-lOafx on 2011 September 
18, 2011 September 21, and 2011 December 7 (phases 
+ 153 d, +154 d, and +186 d), respectively (Tabled]). 

The images were flat fielded, sky subtracted, and 
stacked using standard tasks in IRAF, including the 
gemini package for the NIRI data, except for the HAWK- 
I data, which were processed using the instrument 
pipeline. The images from each instrument were then 
calibrated using the same 2MASS stars in each field, ex- 
cept for the Y image, which was calibrated using archival 
NIRI zeropoints. The H and Ks host galaxy fluxes were 
then subtracted numerically from each datapoint in those 
filters. For the Y data point, we interpolated the host 
galaxy flux between the measured values in j/pi and J to 
subtract off the host galaxy contribution (only a small 
correction). We were able to perform image subtraction 
on most of the J images using the ISIS software package 
with the late-time FourStar image as a template, but 
the others had the host J flux subtracted numerically. 
The final host-corrected photometry points are shown in 
Figure 13] and listed in Tabled] including the calibration 
uncertainty in the errors. 

2.3. Spectroscopy 

After noticing the unusually red rpi — ipi colors of PSl- 
lOafx, we immediately triggered optical spectroscopy us- 
ing GMOS-N. A pair of dithered 1800 s observations were 
taken on 2010 September 11.26 (phase —2.8 d) with the 
R400 grating and OG515 order-blocking filter to cover 
the wavelength range of 5400-9650 A. We used stan- 
dard tasks in IRAF to perform basic two-dimensional 
image processing and spectral extraction. We used our 
own IDL procedures to apply a flux calibration and cor- 
rect for telluric absorption based on archival observations 
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of spectrophotometric standard stars. The observations 
were performed at a mean airmass of 1.5, but the 1"- 
wide l ong slit w as oriented within 20° of the parallactic 
angle (jFilippen ko 1982), so the overall spectral shape is 
reliable at these red wavelengths. 

A second epoch of spectroscopy was obtained on 2010 
September 26.04 (phase +3.4 d) using GMOS-S. The ob- 
servations and instrumental setup were similar to the 
first epoch, except that a redder grating tilt was used 
to cover the range 5880-10100 A. This spectr um has a 
lower signal-to- noise ratio (S/N) than the first one, but 
the main SN features are still present. We attempted a 
final epoch of spectroscopy using GMOS-S in nod-and- 
shuffle mode on 2010 November 5, but PSl-lOafx had 
faded significantly, so we only detected a faint contin- 
uum with a hint of SN features. We do not consider this 
spectrum further. 

We observed PSl-lOaf x with the F olded-port In- 
fraRed Echellette (FIRE; iSimcoe et all [2008) spectro- 
graph on the Magellan Baade telescope on 2010 Septem- 
ber 18.17 (phase -fO.l d), using the high throughput, low- 
resolution prism mode to cover the range of 0.8—2.5 /im. 
The spectral resolution is a strongly decreasing function 
of wavelength, but is i? ~ 500 in J. Six 150 s dithered 
exposures were obtained for a total of 900 s of on-source 
time. Despite the short exposure time, the continuum 
of PSl-lOafx is well detected blueward of ~1.6 /im. We 
reduced and combined the spectra using the FIREHOSE 
package, including a correction for telluric absorption ob- 
tained from observa tions of the AOV sta r HD 208368 us- 
ing the algorithm of iVacca et al.l (|2003l ). 

We obtained a MMIRS spectrum on 2010 September 
28 (phase -1-4.2 d), using the zJ filter and J grisni to 
cover the range 1.0—1.34 /im. A total of 90 minutes were 
spent on source in a series of dithered 300 s exposures. 
We used standard tasks in IRAF to combine and extract 
the spectra. 

The final optical and NIR spectra are shown in Fig- 
ure m All of our optical spectra exhibit an emission line 
near 8903 A, which we identify as [O II] A3727 emission 
from the host galaxy, as well as absorption from the Mg II 
A2800 doublet at the same redshift of z = 1.388, which 
we adopt as the SN redshift. 

2.4. Radio Observations 

Popular models for SLSNe do not predict de- 
tectable radio emissio n at high redshift, as discussed by 
iChomiuk et al.1 ()201in . However, given its unusual prop- 
erties, we observed PSl-lOafx on 2010 September 12.2 
(pha se —2.4 d) with th e Karl G. Jansky Very Large Ar- 
ray ( Perlev et al.ll20lTI ) as part of our NRAO Key Science 
Project "Exotic Explosions, Eruptions, and Disruptions; 
A New Transient Phase-Space," with 256 MHz of band- 
width centered at 4.96 GHz. The data were reduced 
with standard task s in the Astron omical Image Process- 
ing System (AIPS: 'Greisen!l2003f). using J2212-f0152 as 
the gain calibrator and 3C48 as the absolute flux density 
calibrator. These observations yielded a non-detection of 
19±18 /fJy. 

3. HOST GALAXY OBSERVATIONS 

We measure a host redshift of z = 1.3883 ± 0.0001 
from fits to the [O II] line in our highest S/N spec- 
trum (phase —2.8 d). A double- Gaussian fit to the 



Table 2 

PSl-lOafx Host Galaxy Photometry 



Filter Observed/Rest-frame Magnitude^ Error Instrument 





Wavelength (A) 


(AB) 






g' 


4825/2020 


24.01 


0.10 


LDSS3 


r' 


6170/2580 


23.69 


0.10 


LDSS3 


«Pi 


7520/3150 


23.43 


0.13 


PSl 


2P1 


8660/3625 


22.75 


0.10 


PSl 


ypi 


9620/4030 


22.29 


0.28 


PSl 


J 


12500/5230 


22.02 


0.14 


FourStar 


H 


16500/6910 


21.93 


0.22 


FourStar 


Ks 


21490/9000 


21.53 


0.25 


FourStar 



^ Corrected for Galactic extinction. 



Mg II A2800 doublet absorption in the same spectrum 
(inset in Figure S]) is blueshifted by 120±12 km s^^ 
relative to the [O II] rest frame. Such an offset is 
typical of rest- fram e UV selected star- forming galax- 
ies at this redshift ()Erb et al.ll2012|) and was also seen 
in ob servations of the SLSN PSl-llbam (jBerger et al.l 
120121 ). This effect has been interpreted to be caused 
by absorption occurring in galactic-scale outflows driven 
by star formation. The rest-frame equivalent widths 
{Wr) of the two lines are iy^(A2796)=1.8±0.2 A and 
Wr(A2803)= 1.6±0.2 A. These v alues are larger than for 
PSl-llbam (jBerger et al.l 1201 2f ). but slightly lower than 
the median of the intrinsic absorbers from g amma-ray 
burst (GRB) host galaxies (|Fvnbo et al.ll2009f) . In addi- 
tion, Mg I A2852 absorption is present, but the UV spec- 
tral slope of PSl-lOafx is so red that the S/N rapidly 
decreases to the blue, making it hard to confirm the 
presence of other expected strong interstellar absorption 
lines, such as Fe II A2600. 

We obtained photometry of the host galaxy of PSl- 
lOafx in eight filters from g' to Kg (Table [2|) to probe 
the host stellar population from the rest-frame NUV 
(~2000 A) to the NIR (-0.9 /zm). The host is well de- 
tected in our deep ipizpiypi PSl template images, which 
we supplemented with the observations from other facil- 
ities described above. Host photometry was performed 
with a consistent aperture of radius 1.7" in all filters. 

W e fit a suit e of sin gle stellar population age models 
from iMarastonI (|2005f ) to the data, assuming a red hori- 
zontal branch morphology, a Salpeter initial mass func- 
tion, and a metallicity oi Z — 0.5 Zq. The best-fit model 
is shown in Figure[5]and has an age of 10^ yr and requires 
a small amount of internal extinction (Ay = 0.4 mag). 
The derived stehar mass is ~1.8 x 10^° M©. The g' 
magnitude of the host galaxy corresponds to a NUV 
continuum luminosity of 1.3x10^^ erg s~^ (after correc- 
tion for internal extinction), which implies a star for- 
mation rate_(SFR)of~ 18 M© yr~^ using the calibra- 
tion of iKennicuttI (|1998f ). We also estimate a consistent 
value of ~13 M© yr"-*^ from the observed [O II ] flux of 
(4.7±0.5)xl0-i^ergcm-2 s"! (|Kennicuttlll99l . These 
SFRs would be a factor of ~2 lower without any correc- 
tion for internal extinction. 

This galaxy is significa ntly more massive than the 
previous hosts of SLSNe. iNeill et al.l (|201lD examined 
the host galaxies of a sample of luminous SNe and 
found them to have generally low stellar masses and 
high specific SFRs. In particular, all of the SLSNe-I 
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Figure 4. Optical (black) and NIR (blue) spectra of PSl-lOafx. The spectra are labeled with their phase in rest-frame days. The O II 
A3727 emission and Mg II A2800 absorption from the host galaxy are labeled, as is the broad Ca II H+K P-Cygni feature from the SN 
itself. The zeropoint of the flux scale is appropriate for the day -2.8/-I-0.1 GMOS/FIRE spectra. The day -I-3.4/-I-4.2 GMOS/MMIRS 
spectra have been shifted down by 2.5 flux units. The gray bars labeled with a ffi symbol represent the regions of strong telluric absorption 
between Y I J and J /H. The inset shows the fit to Mg II A2800 absorption doublet. 
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Figure 5. Host galaxy SED. The black data points represent the 
observed host galaxy photometry (Table [2]l , corrected for the de- 
rived Ay = 0.4 mag of internal reddening, and plotted versus ob- 
served wavelength. The filled red squares show the expected fluxes 
of the best-fit model (red line) in each of the observed bandpasses. 
The inset shows the 1, 2, and Scr contours for the best fit in the 
host reddening-age plane. 



in their sample had dwarf hosts, with a median stellar 
mass of '--^2x10* M©. By comparison, PSl-lOafx has 
a significantly more massive host with a lower specific 
SFR (~8xl0~^° yr"^) than their sample (with a me- 
dian of 3xl0~^ yr~^). Only the luminous host galaxy 
of the SLSN-II 2006gy exceeded our estimate for the 
stellar mass of the host of PSl-lOafx. Objects dis- 
covered subsequently have continued this trend, having 



either dwarf host galaxies or no host at all detected 
to date despite deep observations (iQuimby et all 120111: 
Chomiuk et all 120111: ILeloudas et al.l l2012t iChen et al.l 
2013; Lun nan et ani2013D . 

While it has been argued that the association of SLSNe 
with dw arf galaxy hosts is driven by their low metallic- 
ity (e.g.. lNeill et al.ll2011| ). the massive host of PSl-lOafx 
appears to be inconsistent with a low metal l icity. For 
example, the relationships of iMannucci et all ()2010[ ) im- 
ply a super-solar metallicity of 12 + log(0/H) = 8.85 for 
our derived stellar mass and SFR, although we caution 
that we cannot measure the metallicity directly with the 
available data. An additional point of contrast is pro- 
vided by the sample of long-duration GRBs, whose oc- 
currence is also widely believed to be associated wit h low- 
metallicity environments (e.g.. iStanek et aT]|2006D . The 
host galaxy of PSl-lOafx is brighter in Kg than any of the 
GRB host galaxies at similar redshifts in the optically- 
unbiased sample of iHiorth et all (1201 2[). although a few 
of the "dark" GRB hosts studied by IPerlev etatl (pOll ) 
are similarly luminous. 

The absolute magnitudes of PSl-lOafx's host prior to 
any internal extinction correction, M ° « —20.3 mag 

•' ' 2800A ° 

and Mb ~ —21.7 mag, are ~0.4 mag brigh t er tha n M*, 
the characteristic magnitudes in ISchechteil (jl976D func- 
tion fits to the field galaxy luminosity functions at this 
redshift (|Gabasch et all [200l iFaber eFall [2001 . The 
Ks magnitude of the host is also near the median of the 
distribution of field g alaxies at this redshif t weighted by 
star formation rate (jKaiisawa eFaLllMol ). Combined, 
these points suggest that the environment of PSI-lOafx 
is representative of typical star-forming galaxies at this 
redshift, while being distinct from the other hydrogen- 
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poor SLSNe. 

The host is marginally resolved in our template im- 
ages, with a full width at half-maximum (FWHM) of 
~1'.'2, or 10 kpc at this redshift. Astrometric alignment 
of difference images taken near maximum light with the 
templates shows that PSl-lOafx is aligned with the cen- 
troid of the host to within O'.'l (0.8 kpc) The other SLSN 
with a massive ho st galaxy, SN 2006gy ()Ofek et al.ll2007l: 
iSmith et al.ll2007t ). is also close to the nucleus, perhaps 
implying an unusual star formation environment. 

4. SPECTROSCOPIC COMPARISONS 

Despite its extraordinary luminosity (see Section 5), 
PSl-lOafx has spectra that more closely resemble those 
of a normal SN Ic than any known SLSN. We show our 
highest S/N spectrum (phase —2.8 d) in Figures [6] and 
[7] along with several comparisons drawn from the liter- 
ature. Relatively few UV spectra of core-collapse SNe 
exist, but the available data allow us to sample a wide 
variety of phenomena. The most prominent feature in 
the PSl-lOafx spectrum is a broad P-Cygni feature with 
an absorption minimum near 3730 A, which we identify 
as the typical Ca II H&K absorption seen in most types of 
SNe, blueshifted by about 16, 000 km s^^. There are sev- 
eral other weaker features in the spectrum in the range 
3000-3500 A. The +3A d CMOS spectrum is noisier, 
but similar, with some evidence that the broad Ca II 
absorption minimum decreased in velocity, although the 
overlapping [O II] emission and strong night sky residuals 
make this uncertain. 

These spectra are in strong contrast to the SN 2005ap- 
like class of SLSNe (|Ouimbv et al.lJMll). PSl-lOky is a 
high-redshift example of this class (^ Chomiuk et al.ir201lD 
and it has a much bluer NUV continuum with no Ca II 
feature (Figure [6]). A reasonable concern is that PSl- 
lOafx only appears to have redder colors due to host- 
galaxy extinction. However, matching the UV spectral 
slope of PSl-lOafx to PSl-lOky requires E{B - V) ^ 
1 mag of extinction (gray line in Figure Et, as suming a 
Galactic reddening curve (ICardelh et al.l1l989l) . This is 
unsatisfactory for two reasons. First, PSl-lOafx is al- 
ready more luminous in the NUV than the SN 2005ap- 
like objects (see below), so A„ « 4.7 mag would imply 
an upward correction of two orders of magnitude to an 
already extreme peak luminosity. Second, the strongest 
UV spectral features of the SN 2005ap-like objects are 
the trio o f features blueward of 3000 A first seen in 
SCP06F6 (jBarbarv et al.ll2009[ ). with the continuum red- 
ward of that being la rgely smooth except for O II features 
(jOuimbv et all 1201 1[ ). PSl-lOafx has fundamentally dif- 
ferent spectral features with no correspondence in the 
PSl-lOky spectrum, so the simplest explanation is that 
it has an intrinsically cooler photosphere. 

The SN Iln 1998S is shown in Figure IHl as a well- 
studied example of a SN undergoing strong CSM in- 
teraction. SN 1998S at this phase showed no broad P- 
Cygni features in the UV, but had a blue continuum with 
narrow absor ption lines froni inter stellar and circumstel- 
lar material (jFransson eFall 120051 ). The SLSNe 2006gy 
and 2006tf are higher-luminosity versions of SNe dom- 
inate d by hydrogen-rich CSM intera ction (|Ofek et al.l 
[2007t ISmith et alJ [2007L [20081 IMol) . These objects 
have redder spectra than SN 1998S, but still exhibit 
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2000 2500 3000 3500 4000 
Rest Wavelength (A) 

Figure 6. SLSNe and other NUV spectral comparisons for PSl- 
lOafx (blue; —2.8 d). None of these objects is a good match to 
PSl-lOafx. The spectrum in gray is the same one of PSl-lOafx 
after correction for E{B — V) = 1.0 mag of reddening (see text 
for details). The other spectra and their phases relative to maxi- 
mum light are: SN 2007bi (+54 d; Gal- Yam ct al. 2009 ), SN 2008e s 
(-1-21 d: Miller et al. 2009), SN 1998S (+14 d; Fran sson et al.l2005l 1, 
PSl-lOky (-2 d;E;homiuk et al. 2011), SN 2002ap (+4 d; data re- 
trieved f rom the Hubble Sp ace Telescope [HSJ] archiv e), SN 2006tf 
(+64 d; ISmith et ani20OT l. and SN 2006gy (+23 d; ' Smith et afl 
I20igi . All literature data are corrected for Galactic and host red- 
dening, if known. 

Balmer lines and deep narrow absorptions that have no 
analog in the PSl-lOafx spectrum. The SLSN 2008es 
had almost feature l ess blue spectra n ear maximum light 
(jMiller et al.l 120091 : iGezari et al1 f2009). but the plotted 
spectrum is from a later epoch when Balmer lines were 
beginning to develop strength in the rest-frame optical 
spectrum. 

The NUV spectrum of the broad-hned SN Ic 2002ap is 
included because the engine-driven SNe associated with 
GRBs have similar s pectra, even if 2 002ap itself was not 
unusually luminous (jMazzah et al.| [2002). SN 2002ap 
has much broader and more blended features than PSl- 
lOafx. Thus, none of these objects in Figure [S] resem- 
bles PSl-lOafx. The SLSN 2007bi has been identified 
as a potential pair-instability SN (jCal-Yam et al.. ,20091 : 
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Figure 7. Normal SNe NUV spectral comparisons for PSl-lOafx 
(blue; —2.8 d). These objects are much better matches than the 
ones shown in Figure|6] The other spectra and their phases rel ative 
to maximum light are: SN 1993J (-1 d; IJefferv et allll99l) . SN 
19941 (+10 d and +35 d; data retrieved from the HST archive ) 
and two SNe la, SN 2011 fe (+0 d; 'Maguire et al.' 2012; Folev'Mig) 
and SN 2011iv (+0.6 d; IFolev et al. 2012). The SNe la and SN 
1993J spectra have been blueshifted to approximately match the 
absorption minima of the Ca II feature. All literature data are 
corrected for Galactic and host reddening, if known. 



lYoung et al.ll2010l ). It is more similar to PSl-lOafx, with 
a definite Ca II feature, but is still significantly bluer 
and has several additional spectral features that do not 
match. 

Instead, the set of spectra from normal SNe plotted 
in Figure [7] provide a much better basis for compari- 
son. The SNe la and SN 1993J (Type lib) spectra have 
been blueshifted by the indicated amounts to approxi- 
mately match the Ca II absorption minima to PSl-lOafx. 
Over the range 2800—4000 A, the spectra in Figure [7] 
are all very similar, with the dominant features being 
due to Ca II in all objects. A notch near 3950 A IS 
from Si II A4130. An emission peak at 2900 A also ap- 
pears to be common to these objects, along with a bluer 
peak near 3150 A. Shortward of ~2800 A, the SNe la 
diverge from the core-collapse objects. This dropoff in 
flux, along with a pair of deep absorption features near 
2450 A and 260 A, is caused by ir on-group elements , 
primarily Co II ()Kirshner et al.lll993HSauer et al.ll2008[ ). 
The lack of these features in PSl-lOafx is evidence that 
the abundance of newly-synthesized ^^Ni (which decays 



to Co) near the photosphere is much lower than in a 
SN la. 

One purpose of these comparisons is to demonstrate 
the apparent broad similarity of the NUV spectra of PSl- 
lOafx to those of SNe coming from three very different 
progenitors: a bare core of a mass i ve st ar (SN 19941: 
iWheeler et "all 119941: llwamoto et"aLl I1994D . a partially- 
stripped massive star in a binary system that retained 
part o f its hydrogen envelope (SN 1993J: Nomot o et al.l 
[l99l iPodsiadlowski et al.l [1991 Imoslev et al.l I1994D . 
and the thermonu clear explosions of white dwarfs (SN e 
2011fe and 2011iv: iNugent et al.l201lHFolev et al.l2012} ). 
SNe 20 life and 201 liv have some spectral differences 
from each other near 3000 A, exhibiting the potential 
for spectroscopic variation in objects with similar pro- 
genitors. Other than the lack of cobalt near the photo- 
sphere, the remaining differences between PSl-lOafx and 
these other objects are relatively minor. 

In addition, the overall UV spectral slope of PSl-lOafx 
is very similar to those of SNe 1993J and 19941, after cor- 
rection of both by E(B — V) = 0.30 ma g of reddening 
()Richmond et al.l (l994': 'Sau er et al.| [200g). The fact that 
the only SNe in Figures [5] and [7] that have similar spectral 
features to PSl-lOafx also have similar SED shapes pro- 
vides further evidence that the extinction of PSl-lOafx 
is not large. If the PSl-lOafx spectra were corrected for 
a large amount of reddening, they would be bluer and we 
would expect different spectral features to be present at 
the higher implied temperatures. 

We further explore the spectra using the SN spec- 
trurn synthesis code SYN++ with SYNAPPS (Th omas et al.l 
[20T1 to fit the near-maximum-light CMOS/FIRE spec- 
tra. These codes are b ased on the sa.me as sumptions as 
the original SYNDW (e.g.. [Branch et al.|[2002t ) treatment of 
resonant-scattering lines in the Sobolev approximation 
above a sharp photosphere that emits as a blackbody 
(BB). The atomic and ionic species are assumed to have 
levels populated in local thermodynamic equilibrium at 
some excitation temperature. SYNDW is a tool designed 
for line identifications and is not a full self-consistent 
spectral model. 

Our best SYNDW model is shown in Figure [8] and con- 
tains seven ions. We assumed equal BB continuum and 
excitation temperatures of 10, 200 K and a photospheric 
velocity of 11,000 km s~^. The deep absorption trough 
near 3700 A is dominated by Ca II H&K. Fe II and Ti II 
combine to produce several of the wiggles observed in 
the range 3000 - 3500 A. Mg II naturally explains the 
observed absorption near 4250 A, while also producing 
absorption shortward of ~2700 A. Si II contributes to 
both the blue wing of the deep Ca II feature and the 
notch near 3950 A. The feature near 3950 A is present in 
the overlap between the CMOS and FIRE data and has a 
consistent strength that is significantly larger than in the 
synthetic spectrum. With SYNDW, we have to be careful 
not to overproduce absorption by the Si II A6355 doublet 
because no strong feature is seen in the data, although 
the S/N and the spectral resolution become very poor in 
the observed-frame H. We introduce Fe III to help re- 
produce the wiggle near 5000 A, but the identification is 
not unique and its presence is not required. In addition, 
we include Cr II to help suppress the UV flux. That ion 
does not contribute sufficiently strongly to any specific 
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Figure 8. SYNOW fit for PSl-lOafx (red) compared to tiie 
— 2.8/+0.1 d GMOS/FIRE combined spectrum. Ions dominating 
major spectral features in the SYNOW model are indicated. The 
gray bars labeled with a ffi symbol represent the regions of strong 
telluric absorption. 

feature for us to make a positive identification, but it 
has been included i n previous SYNOW models of SNe Ic 
([Millard et al.lll999l ) because of its expected strength in 
either helium or carbon/oxygen- r ich SN ejecta at these 
temperatures (|Hatano et al ill99"9l ). 

We also consider several other species that were not in- 
cluded in the final fit. As discussed above, Co II produces 
strong NUV absorptions in SNe la that are not present in 
PSl-lOafx. Na I D is commonly seen in all types of SNe, 
but at this redshift it would fall in the strong telluric 
absorption between J and 7J, so we have no constraint 
on its presence from the observations. We also tested 
Sc II b ecause of its use in prior SYNOW analyses of SNe lb 
([Branc h et al. 2002). Scandium results in an additional 
absorption minimum near ^3500 A similar to one seen 
in our data, but does not produce a sufficiently strong 
improvement in the model to justify its inclusion, partic- 
ularly given the expectation that it should only become 
impo rtant at lower ejecta temperatures (jHatano ct all 
Il999| ) . There is also an apparent absorption feature near 
5300 A that is not fit by the model. S II can produce 
a feature near that wavelength, but the identification 
is not certain. Lines of hydrogen or helium would be 
very important if present, but there is no evidence for 
either, although the low S/N of our NIR data precludes 
strong statements about helium. The SN 2005ap-like 
SLSNe have lines of singly-ionized carbon and oxygen i n 
their spectra at early times (e.g., ' Quimbv et al.|[2"011[ ). 
with possible contributions from doubly-ionized CNO el- 
ements ([Pastorello et al.ll201Q ). but those lines do not 
appear to be present in our PSl-lOafx spectra. 

Overall, we can reproduce the shape of the spectrum 
in the rest-frame optical reasonably well, but the NUV is 
more problematic. Our synthetic spectrum has features 
at many of the same wavelengths as the real data, but 
the overall shape in the NUV is less well fit. Adding 
additional ions to the fit might potentially help the 
UV shape of the synthetic spectrum, but would not be 
well motivated or add any additional insight. Alter- 
natively, the radiative transfer in the UV is known to 
be complex and we could be encountering the limita- 



tions of SYNOW. The most important point is that the 
ions we have used are typical for generic fits to SNe 
la and Ic spectra (e.g., [Branch et al.ll200^ . but differ 
from those id entified in most SLSNe-I spectra near max- 
imum light (Quimbv et al.' 2011; Chomi uk et al.l 1201 ll : 
iLeloudas et al. 2012: Lunnan et a l. 2013). 

5. LIGHT CURVE AND ENERGETICS 

Existing light curves of SLSNe are heterogeneous, 
and the objects have been found at a wide range of 
redshifts, making d irect comparisons d i fficult. Follow- 
ing previous work (iQuinibv et al.l l 201ll : iChomiuk et al.l 
[201lt ILunnan et al.l 120131) . in Figure [H] we compare the 
rest-frame u light curve of PSl-lOafx (corresponding 
to observed- frame z-pi) to several of the most lumi- 
nous SN 2005ap-like SNe in the literature. In addi- 
tion, we include the hydrogen-rich SN 2008es because 
its light curve is similar to the others and it did not 
devel op strong Balmer l ines until well afte r maximum 
light ([Miller et al.l [20091: iGezari et al.] [20091) . For com- 
parison, we also plot the U light curves of the broad- 
lined SNe Ic 1998bw and 2002ap, two objects of nor- 
mal luminosities (Galama ct al. 1998; Foley ct al. 2003). 
Given the uncertain SEDs of many of these objects, we 
do not perform a full fc-correction, but instead correct 
the observed AB magnitudes in filters whose effective 
wavelengths are close to u in the observed frame by 
m — M = 5 log(dL(z)/10 pc) — 2.5 \og{l-\-z), where d^ is 
the luminosity distance. To be conservative, we assume 
that PSl-lOafx is unaffected by reddening here and in all 
following discussion, despite the strong Mg II absorption 
along the line of sight and the inference that the host 
galaxy SED requires some degree of internal extinction. 
We presented some evidence above that the reddening 
to the SN is not large, but it is possible that we have 
underestimated the peak luminosity of this object. 

It is apparent from Figure [S] that at peak PSl-lOafx 
was slightly more luminous in u (M„ k, —22.3 mag) 
than any previous SLSNe. However, by ^10 d on either 
side of maximum light, PSl-lOafx falls below all the oth- 
ers. This emphasizes the high peak luminosity and fast 
evolution of the light curve. These objects have rather 
different SEDs, so comparisons of the energetics require 
constructing a bolometric light curve. 

To proceed further, we first construct SEDs at three 
epochs near the times of our NIR observations, at phases 
of -1-3, -f 12, and -f 24 d. We fit third-order polynomials to 
the observed photometry in each rpiipizpi filter around 
each SED epoch to estimate the flux. We then estimate 
the error bars on these fitted fluxes by generating Monte 
Carlo realizations of the data where we repeatedly ran- 
domly adjust the observed fluxes by drawing from a nor- 
mal distribution having a width of the measured errors, 
refitting, and determining the variance of the interpo- 
lated flux at the desired epoch. We perform a similar 
procedure whenever we need to interpolate fluxes to a 
common epoch. For the SEDs at phases of -1-3 and -1-24 d, 
we use the actual measured ypi and NIR fluxes. For the 
-1-12 d SED, the observed NIR fluxes were quite noisy, 
so we again interpolate the fluxes in J and H by fitting 
second-order polynomials to the observed light curves to 
generate the points on the SED. The resulting SED evo- 
lution is shown in Figure 1101 

We then fit BB spectra to the SEDs at each epoch. 
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Figure 9. Absolute u light curves of PSl-lOafx and other SNe. 
The observed bands and wavelengths in the rest frame are: PSl- 
lOafx (zpi, 3525 A), SCP 06F6 (F850LP, 3883 A: IBarbarv et al.1 
[2009D , PSl-lOky (ipi, 3848 A: IChomiuk et al.l 120111) . PSl-lOawh 
(ipi, 3944 A: IChomiuk et al.1 120111) . PTF 09cnd (rest-frame u 
from'Quimbv et al."2011D, SN 2008es (B/g, 3604 A: IMiller et"an 
[2009; Gezari et al. 2001), SN 2010gx/PTF lOcwr {g, 3770 A; 
[Fastorello et al. 2010; Ouimbv et al.l2011l) . SN 1998bw (U, 3663 A; 
[Galama et al. 1998), and SN 2002ap (U, 3663 A; IFolev et al.l2003l) . 

Our -1-3 d SED has the most complete information. The 
dashed hne in FigurefTOlshows the result of a BB fit to the 
full SED. The derived BB temperature is Tbb ~ 6800 K, 
but clearly the curve is a poor fit to the data at rest 
wavelengths longward of ~5000 A. If instead we restrict 
the fit to j/pi and the NIR bands, the BB temperature is 
not well constrained (Tbb = 12, 000 ± 2, 000 K) due to a 
lack of points blueward of the peak. The formal best fit 
is plotted as the solid blue line in Figure [TU] and provides 
a much better fit to the NIR data points (by construc- 
tion), but significantly overestimates the rest-frame UV 
flux. This is to be expected if line blanketing in the NUV 
provided by iron-peak elements supresses the UV flux, as 
is typical for SNe Ic. For example, the shape of the SED 
of the normal SN Ic 19941 near maximum light (Fig- 
ure (TUl) is simila r to that o f PSl -lOafx, but the detailed 
model fits by Sauer et all (|2006D find an underlying BB 
temperature of 9230 K at the same epoch. 

The data at phases -1-12 and -|-24 d are less complete, 
so we fit a simple BB to all of the points on the SED 
and do not correct for line blanketing, which would only 
be expected to increase as the ejecta cool. The derived 
BBs show the ejecta clearly cooling from '^12, 000 K at a 
phase of -^3 d to '-5500 K at -^24 d, ahhough Tbb should 
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Figure 10. SED of PSl-lOafx at three epochs, phases of +3 (blue 
triangles), -|-12 (green diamonds), and -f24 d (red circles). The 
solid lines are BB fits to the SED at each epoch. At a phase of -|-3 d, 
the blue dashed line is a fit to all data points, while the solid line 
excludes gpi»'pijpizpi from the fit. Reported BB temperatures 
are in the rest frame. The magenta line shows the SED of SN 
19941 near maximum light shifted to the redshift of PSl-lOafx and 
arbitrarily scaled to match the fiux level. The 19941 data poi nts are 
the U BVRI photometry on 1994 April 5.39 from Richmo nd et al.l 
((199^ corrected for E{B - V) = 0.3 mag liSauer^tal.l 12009 1 of 
reddening. 
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Figure 11. Observed color evolution of PSl-lOafx. At this red- 
shift, 2pi is near rest-frame u band and J is near V, while jpi is 
in the NUV. Note the flat ipi — zpi color prior to maximum light. 

only be regarded as lower limits to the effective temper- 
ature (Toff) on the later two dates. Independent of the 
details of line blanketing, we can directly see the SED be- 
coming redder, and hence implying cooler temperatures, 
by looking at the color evolution. In Figure 111! we con- 
struct ipi—zpi and zpi — J color curves. In each case, we 
interpolate the light curve in the bluer band to the epochs 
of observation for the noisier redder band. The ipi— zpi 
color is consistent with a constant (0.30±0.05 mag) be- 
fore maximum light, but both colors became redder with 
time after maximum. 

We use the BB fits to construct a bolometric light 
curve, starting with the near-maximum-light SED (phase 
+3 d) shown in Figure [TOl A trapezoidal integra- 
tion of the observed fluxes from rpi to K gives a 
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minimum bolometric luminosity (Lboi) near peak of 
(3.6±0.2)xl0*'* erg s~^. To account for flux emitted 
redward of our observations, we add the tail of the BB 
and extend the trapezoidal integration to long wave- 
lengths, which provides only a 14% upward correction 
to the total flux. We repeat the process for the other 
two SED epochs. After including the BB tails, the de- 
rived bolometric luminosities are (4.1±0.2), (2.2±0.2), 
and (l.l±0.3)xl04'' erg s"^ at phases of +3, +12, and 
-1-24 d, respectively. 

Although we lack NIR data at other epochs, we can use 
these fits to estimate the bolometric light curve in combi- 
nation with the well-sampled PSl light curves. For each 
of our three SED epochs, we define a multiplicative bolo- 
metric correction factor to obtain Lboi from vL^j in the 
zpi band. Motivated by the lack of strong color evolu- 
tion before maximum light, we apply the bolometric cor- 
rection factor from -|-3 d to all pre-maximum zpi data 
points. For the post-maximum data, we use a smoothly 
varying time-dependent correction factor to evolve from 
the -1-3 d value to the -1-24 d value. We repeat the process 
for the ipi band, except that we use a smaller (but con- 
stant) bolometric correction factor before maximum light 
than the derived -1-3 d value to better match the zpi re- 
sults and account for the small amount of i-pi—z-pi color 
evolution between -f3 d and the pre-maximum data. 
The combined bolometric light curve is plotted in Fig- 
ure [12] A third-order polynomial fit to the peak of the 
bolometric luminosity curve gives a maximum value of 
(4.1±0.1)xlO'*'* erg s'S equivalent to 1.1 xlO" Lq, not 
including the uncertainty in the bolometric correction 
factors (or potential extinction). 

To check that our derived maximum bolometric lumi- 
nosity of PSl-lOafx is realistic and that our procedure 
for extrapolating the SED to the NIR has not intro- 
duced any significant error, we can take advantage of 
the fact that the shapes of the spectra we have of PSl- 
lOafx and the derived SEDs are very similar to those 
of a SN la near maximum light (cf. Figure [7]). One 
of the most distant spectroscopically confirmed SNe la 
known, HST04Sas, is at z = 1.39 (Riess et al. 2007), con- 
sistent with the redshift of PSl-lOafx. Interpolating the 
F850LP (-Z band) light curve of HST04Sas to a phase 
of zero days produces an estimate of the observed mag- 
nitude of mF85QLp{t = 0) = 25.5 mag (AB), while fits 
to PSl-lOafx show that zpi{t = 0) = 21.7 mag. Assum- 
ing that the similarity of the SEDs requires no further 
correction factors and a typical SN la peak bolometric 
lum inosity of ~1.2xl0''^ erg s~^ (e.g., iStritzinger et al.l 
12006) ■ the 3.8 mag difference between the objects corre- 
sponds to a peak bolometric luminosity for PSl-lOafx of 
Lho\ ~ 4x10^** erg s~^, in excellent agreement with our 
value derived above. 

In fact, the peak of Afboi = —22.8 mag is more lu- 
mi nous than any o f the hydrogen-poor objects collected 
by iGal-Yaniil (|2012rF^ . This statement is somewhat un- 
certain due to the heterogeneous and fragmentary na- 
ture of the available data for SEDs and bolometric cor- 
rection factors for the ob j ects in the l i teratu re. SNe 
2n08es (IMiller et al.l [2nn9l: IGezari et al.l l2009h. 2005aD 
(|Quimbv et al.l I2007D . and SCP06F6 (jBarbarv etal] 

For consistency, we note that using the same cosmology as 
Gal- Yam would result in Mjjqj = —22.9 mag. 
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Figure 12. Bolometric light curve of PSl-lOafx. The observed 
ipi (circles) and zpi (squares) light curve points were multiplied by 
a time-dependent bolometric correction factor (see text for details). 
The plotted error bars only represent the nominal photometric er- 
rors and do not include any contribution from the systematics of 
converting to bolometric luminosities. The solid blue line shows a 
fit to the early-time flux that rises as . The red dashed line shows 
a model for a SN powered by the dipole spindown energy from a 
magnetar, following Kasen & Bildsten (2010). We assume an opac- 
ity of 0.2 cm^ g~"^, 1 Mq of ejecta, 10^-^ erg of kinetic energy in 
the initial SN, an initial magnetic field strength of 6x10^'* G, and 
an initial magnetar spin period of 2 ms. This model predicts high 
temperatures and velocities that are in contradiction of the data 
despite the apparent good fit. See Section 6.2 for further discus- 
sion. 

I2009t iQuimbv et al.|[20Tll ) all have similar peak bolomet- 
ric luminosities to PSl-lOafx, but the exact comparison 
depends on the treatment of the bolometric correction 
for these UV-luminous sources. The only SLSN of any 
type with a clearly higher peak lumin osity is the peculia r 
hydrogen-rich transient CSS100217 (|Drake et all 120 111) . 
However, the SN nature of that object is still question- 
able due to possible confusion with the active galactic 
nucleus of its host galaxy. 

In addition to the high peak luminosity of PSl-lOafx, 
another striking fact about the light curve is the rapid 
time evolution. Despite the rolling survey nature of PSl, 
our first detection is at a phase of only —7.4 d, with 
non-detections prior to that (Figure [T2|) . It is tradi- 
tional in SN studies to find the rise time by fitting a 
fireball model of the form iboi oc for t > to a.t early 
times, where tg is the explosion time. For PSl-lOafx, 
this gives a. to = -11.8 ± 0.7 d (Figure [HI). Allow- 
ing the power-law index to vary gives a best fit with 
Lboi oc t0-9±0-5 and = -8.7 ± 1.9 d. Regardless 
of the exact parameterization, the rise in zpi is a fac- 
tor of 3.5 in fiux in the 7.4 d before peak. This ap- 
parent rise time is much faster than the fastest known 
SLSN e (Figure M IQuimbv et al.|[20Tlt iPastorello et al.l 
[Mot iLeloudas et al.ll2012l: iLunnan et al.ll2013[ ). 

The decay timescale is also unusually fast. The FWHM 
of the bolometric light curve is only 18 d. The usual 
parameterization of SN light curve shapes. Amis, is the 
decline in magnitudes in 15 d after maximum light. The 
bolometric light curve of PSl-lOafx has Ami5=0.95 mag. 
By comparison, only SN 19941 has a faster Ami^(R) in 
the co llection of SNe Ib/c light curves of iDrout et al.l 
(|2011f) . although several additional objects had Ami5{V) 
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that are comparable. Figure IHl demonstrates that in u, 
PSl-lOafx also has a faster light curve decay than SLSNe 
of similar peak luminosity. 

Integrating the observed bolometric light curve from 
the first detection at a phase of —7.4 d to the last de- 
tection at +24 d provides a lower limit on the emit- 
ted energy of '^7x10^° erg (over only 31 d). Reason- 
able extrapolations of the light curve to later times will 
increase the total to '^10^^ erg, comparable to most 
other hydrogen-deficient SLSNe, but not notably high 
(e.g.,[Quimby ct al. 2011; Chomiuk ct al. 2011). The fast 
timescale of PSl-lOafx compensates for its high peak lu- 
minosity. 

6. COMPARISON TO MODELS 

The difficult constraints placed by the rise time on 
any model can be seen from the inferred emitting ra- 
dius. Assuming T^s ~ 11,790 K as derived above, the 
BB radius (-Rbb) of the photosphere at maximum light is 
5.5 X 10^^ cm. This is larger than for other hydrogen-poor 
SLSNe, which typically have inferred radi i of ~10^^ cm 
(jOuimbv et all 1201 It IChomiuk et al.ll201lD . as expected 
from the higher luminosity and lower Tbb of PSl-lOafx. 
If Tcff is actually as low as the color temperature of 
~6800 K, then the inferred radius at maximum light in- 
creases to 1.6x10^^ cm. 

The most conservative set of numbers at maximum 
light (i?BB=5.5xlO^^ cm, d) requires an ex- 

pansion velocity of ~50,000 km s~^ if to is the true 
explosion time, in obvious conflict with the observed 
photospheric velocities seen in the spectrum. Material 
moving at the SYNOW-derived velocity of 11,000 km s~^ 
will take 58 d to reach that radius, implying an explo- 
sion ^50 d before our first detection. There is an ob- 
served example of such an effect with the SLSN 2006oz. 
It exhibited a precursor plateau with a luminosity of 
~2xl0^'^ erg s~^ pri or to rising to the m ain peak of 
the light curve, which iLeloudas et al.l (j2012D interpreted 
as being due to a recombination wave in an oxygen- 
rich GSM. A weighted average of the PSl-lOafx non- 
detections between phases —15 and —8 d produces a 
mean luminosity of (0.2±1.2) x 10^'^ erg s~^. At the 2a 
level, this non-detection is consistent with a hypothetical 
plateau being present at the level seen in SN 2006oz, but 
unobserved due to the limiting magnitudes of our survey. 

At phases -1-12 and -1-24 d, the iboi quoted above and 
Tbb shown in Figure [TUl combine to produce BB radii of 
1.4 and 1.3x10^^ cm, respectively, although these are 
likely overestimated due to the underestimate of T^g 
caused by line blanketing. Notably, these values are con- 
sistent with the photosphere begining to recede through 
the ejecta (in both the physical and conioving senses) as 
they become optically thin. By contrast, the BB radii 
inferred for the interacting SLSNe 2006tf and 2006gy in- 
creased linearly with time at the same rate as the veloc- 
ities measured from the spect ra, consistent with an ex - 
panding optically-thick shell (|Smith et al.l I2008L l2010t ). 
Although the measurement errors make it less certain, 
this also appe ars to be true for at least some SN 2005ap - 
like SNe fe.g.. IChomiuk et al.ll20Tl ILunnan et al.ll2013D . 

6.1. Models Powered by Interaction 

The kinetic energy of a powerful SN explosion provides 
an attractive reservoir of energy to power SLSN light 



curves, but only if that energy can be efficiently con- 
verted to radiation. For PSl-lOafx, this already requires 
an initial SN energy of at least a fewxlO^^ erg in order 
for the object to radiate 10^^ erg and still have ejecta ve- 
locities as high as those observed. This generally requires 
a CSM with a mass comparable to that of the ejecta and 
a radius comparable to the photospheric radius at max- 
imum light in order to maximize e fficiency and not lose 
energ y to adiabatic expansion (e.g. JGinzburg &: Balbera 
|2012| ). Another attraction of interaction models is that 
the SN can explode well before the start of the optical rise 
of the SN if the CS M is sufficiently dense to trap radia- 
tion from the shock (jChevalier fc Irwinll201lD . This could 
greatly alleviate the conflict described above between the 
rise time, ejecta velocity, and photospheric radius. 

The required CSM densities in these models are much 
larger than for well-studied examples at lower luminos- 
ity. iChevalier fc IrwinI (j201ir ) have identified two regimes 
of interest when the CSM is this dense. If the radiation 
from the shock can escape while the shock is still interact- 
ing with the CSM, then the light curve will be broad and 
exhibit a slow decay after maximum light. Their model 
works well for the light curves of SLSNe like 2006gy. An 
additional consequence of this model is that there should 
be spectral signatures of the ongoing CSM interaction, 
such as emis sion lines from the unshocked CSM (e.g., 
iChugail 120011 ). Other luminous SNe Iln in the litera- 
ture, such as SNe 200 8fz (Drake et al.ll2010t ) and 2008am 
(jChatzopoulos et al."2011), with their blue continua and 
strong Balmer emission lines, can also be fit by similar 
models. 

The much faster light curve of PSl-lOafx and lack of 
strong emission lines imply a different regime, one where 
the outer radius of the CSM is comparable to or smaller 
than the diffusion length of photons from the shock. In 
this framework, the optical light we detect is then as- 
sociated with the breakout of the shock from the outer 
radius of the CSM. In the context of a dense wind with 
a density profile that scales with radius as p = Dr~^ 
out to a maximum radius of Rw, where Z? is a den- 
sity normalization parameter, the tot al mass of the wind 
is M = AttRwD. In the models of IChevalier fc IrwinI 
(l20Tll) . D can be estimated from the diffusion timescale 
as D = 1.28 X 10^^ (t^/days) g cm~^, if we assume an 
ionized He-rich composition. We can approximate the 
diffusion timescale td by the rise time and set R^ to the 
BB radius at maximum light to produce an estimate of 
the wind mass of 5 — 10 M,-:^. These val ues are only 
approximate, as lGinzburg fc Balberd (|2012l ) have shown 
that in the regime of interest, the analytic approxima- 
tions are outside the range of validity and require true 
hydrodynamic calculations. 

Models with shock breakout occurring in a dense, trun- 
cated wind can produce bolometric light curves with the 
desired peak luminosities and time scales, but it is un- 
clear if this can be done consistently with the other in- 
formation about the object. Ginzburg & Balbcrg (201^ 
were able to approximately match the light curve shapes 
of the SLSNe 2010gx and 2005ap, but the derived Tbb 
were hard to match to the observations. A theoretical 
caveat is that the observed color temperature could de- 
viate from Tpff if the shock is n ot in equilibrium when it 
breaks out (jNakar fc Saril[20Tol ) . This deviation is caused 
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by the nature of the opacity, and usually implies that the 
color temperature is hotter than the true T^g (e.g., see 
IMoriva et al.ll2013l for models of 2006gy) , which would in 
turn require that the photosphere be at even larger radii 
than we have found. This effect might then conflict with 
the fast timescale of the light curve, which otherwise is 
indicative of a more compact C SM (jChevalier fc Irwinl 
120111: iGinzburg fc BalbCT3l2012D . but detailed modeling 
with radiative transfer is needed to correctly evaluate 
these concerns. 

Further, a key component of the ri sing SN luminosity 
in the shoc k breakout models of both | Chevalier fc Irwinl 
((20T1 and IGinzburg fc Balberd (IMa is that the tem- 
perature rises to a maximum near the peak of the light 
curve and declines thereafter. Although we lack full SED 
information on the rise of the light curve for PSl-lOafx, 
the observed ipi — zpi color shows no evolution prior to 
the peak of the light curve fFigure[TT|). Despite the inter- 
pretation uncertainties introduced by the possibilities of 
line blanketing or deviations from equilibrium, the evo- 
lution of Tbb after maximum light results in a noticeable 
reddening of the zpi — zpi color. If Tbb were rising by a 
similar amount while the light curve rose to its peak, we 
would have been able to detect it. 

A separate uncertainty is the origin of this CSM, 
which could not be produced by steady stellar winds 
due to the high required mass-loss rates and the un- 
explained truncation at i?^,. Most theoretical work 
has focuse d on winds, but the niaterial could also be 
in a sh ell (ISmith fc McCrayl 120071 : IMoriva fc Tominagal 
|2012[) . IWooslev et al.l ()2007D invoked a pulsational pair 
instability to produce shells of matter at the appropri- 
ate distances, but these require very massive st ars. Other 
possi bilities include common envelope ejection (jChevalierl 
|2012| ) and instabilities in the l ate stages of stellar evolu- 
tion (|Quataert fc ShioddlMl) . 

Another challenge for CSM interaction being the dom- 
inant contributor to the luminosity of PSl-lOafx is the 
appearance of the spectra. Interaction-powered SNe con- 
vert the kinetic energy of the explosion into UV/optical 
continuum radiation, frequently reprocessing some of the 
continuum into strong emission lines from the dominant 
constituents of the CSM external to the expanding shock. 
The spectral comparisons in Figure [H] show several ex- 
amples of objects consistent with this picture. How- 
ever, PSl-lOafx is even more luminous than those ob- 
jects and yet shows P-Cygni features originating in the 
SN ejecta, including a deep Ca II absorption compara- 
ble to that seen in SN 19941. The lack of strong emis- 
sion lines from the CSM can be explained if the mate- 
rial has a sharp outer radius, but interaction sufficiently 
luminous to dominate the light curve should dilute and 
weaken spectral features from the eje cta due to the "top- 
lighting" effect (jBranch et al.l [2000l ) . if not completely 
hide emission from the unshocked ejecta. For example, 
PTF 09uj is a candidate for being a lower-luminosity ver- 
sion of a shock breakout through a dense (hydrogen-rich) 
wind and, as expected, its spectra arc very blue and lack 
strong P-Cygni features (Ofek et al. 2010). It is not clear 
how strong CSM interaction could dominate the luminos- 
ity of PSl-lOafx without leaving spectral signatures. 

In summary, existing CSM interaction models can 
match the gross photometric properties of PSl-lOafx 
(peak luminosity, timescales) , but only if the models have 



the freedom to assume the necessary CSM structure. The 
origin of such unusual CSM is not explained a priori in 
these models. In addition, the details of the SED, color 
evolution, and emitted spectra appear to be in conflict 
with the observations, although more detailed radiative 
transfer calculations are necessary to be confident in the 
model predictions, which are largely based on BB as- 
sumptions. 

6.2. Internal Energy Sources 

Normal hydrogen-deficient SNe have light curves pow- 
ered by the diffusion of energy deposited by the ra- 
dioactive decay of '^^Ni and ^^Co. The light curves 
of most SLSNe-I appear to fade too rapidly rela- 
tive to their pea k luminosities to have radioactively- 
powered ejecta (IQuimbv et al.l 120111 : iPastorello et al.l 
l2OT0t iChomiuk et al.l 1201 ID . and PSl-lOafx is no excep- 
tion. With a rise time of 12 d, the peak bolometric lu- 
minosity requires a nickel mass of ^^14 Mq. However, 
making st andard assumptions about the structure of 
theejecta (|Arnettlll983: iValenti et aLll2008l: iDrout et al.l 
120111 ) allows us to use the photospheric velocity of 
11,000 km s~^ and light curve timescale of 12 d (assum- 
ing that the diffusion timescale is comparable to the rise 
time) to estimate the ejecta mass as ~2 Mq, clearly in 
contradiction of the high nickel mass. 

Another potential power source suggested for SLSNe 
is the spindown energy from a ne w ly-born magnetar 
iMae da ct al. 2007; Woosl^ 120101: IKasen fc Bildste3 
120101: iDcssart et all I2012af) . The high luminosity and 
fast rise of PSl-lOafx stretch the existing models to the 
bounds of plausib il ity. I nspection of Figures 4 and 5 of 
IKasen fc BildstenI (|2010f l shows that the light curve of 
PSl-lOafx requires a magnetar solution to have a very 
fast initial spin period (p « 1 — 2 ms, near the breakup 
speed of a neutron star), a low ejecta mass (~1 Mq), 
and a high magnetic field {> /ewxlO^'' G). This is 
because the high luminosity requires a high spindown 
power, while the fast light curve requires both fast spin- 
down times and short diffusion timescales. We exam- 
in ed magnetar-powered ligh t curves using the formalism 
of IKasen fc BildstenI (j2010( ) and found parameters that 
produced light curves that approximately matched the 
shape of the bolometric light curve of PSl-lOafx. One 
such model is plotted in Figure [T^ 

However, it is not clear whether this fit should be re- 
garded as more than a numerical curiosity because the 
physical assumptions underlying the model break down. 
In the parameter range appropriate for PSl-lOafx, the 
magnetar dipole spindown energy (5x10^""^ (p/2 ms)~^ 
erg) dominates over the initial SN kinetic energy (as- 
sumed to be 10^1 erg). The IKasen fc BildstenI ()2010l ) 
models make an ansatz that the entire spindown energy 
of the magnetar is simply thermalized spherically at the 
base of the SN ejecta and deposited into the internal 
energy. Because this process accelerates and heats the 
ejecta, the final kinetic energy, the peak luminosity, and 
the temperature are coupled. The model plotted in Fig- 
ure [12] predicts expansion velocities of >25,000 km s~^, 
in contradiction of the observed spectrum. In addition, 
the model predicts a radius at maximum light (from the 
ejecta velocity and rise time) that when combined with 
the luminosity implies a BB temperature in excess of 
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20,000 K (with strong evolution at other times). This 
contradicts the observed Tbb near maximum hght and 
the constant red ipi— zpi color during the rise of the 
light curve. We find that increasing the assumed SN ki- 
netic energy relative to the spindown energy does not 
alleviate these issues because of the difficulty in simul- 
taneously satisfying the constraints from the expansion 
velocity seen at maximum light (which implies a nor- 
mal energy-to-mass ratio, and hence high mass if the ini- 
tial SN energy is high) with the fast diffusion timescale 
(which requires a low mass). 

The models of lWooslevI (|2010f ) start from massive star 
progenitors and proceed through the explosions, but they 
only probe slower rotating, and hence less e nergetic, mag- 
netars than those required for PSl-lOafx. iDessart et al.l 
(j2012a[ ) also examined explosion models for Type Ib/c 
SNe with the addition of a central energy source that 
mimicked the effect of a magnetar. They followed the 
effects of a central energy source on the ejecta produced 
by the explosion models and generated light curves and 
spectra. If the central energy source deposited ~10^^ erg 
of energy, the resulting SNe had high velocities and were 
blue at maximum light, which appears promising for ex- 
plaining the SN 2005ap-like class of SLSNe, but is very 
different from PSl-lOafx. As noted above, the light curve 
of PSl-lOafx requires even more energy input from the 
magnetar. Therefore, we disfavor magnetar models at 
t his time. 

IDexter fc KasenI (j2012) have proposed a scenario that 
shares some similarities with the magnetar model in that 
their model has a post-explosion internal energy source, 
in this case provided by fallback accretion onto a newly 
formed compact remnant, likely a black hole. Their 
suite of models can only match the fast timescale and 
high peak luminosity of PSl-lOafx if the progenitor is 
a blue supergiant and the initial SN explosion energy 
is extremely low with an ejecta mass of ~1 M©. Their 
hydrogen-poor models do not reach the same peak lumi- 
nosities, which is a problem for PSl-lOafx. In addition, 
these models have the same problem as the magnetar 
ones in simultaneously satisfying the high peak luminos- 
ity and relatively cool photospheric temperature. 

6.3. Asphericity 

In all of the preceding discussion, we have necessar- 
ily assumed spherical symmetry in order to convert our 
observations into intrinsic quantities. However, core- 
coUa pse SNe are all known to be aspherical at some level 
(e.g., I Wang fc Wheeleiil2008D . The most basic effect of 
asphericity is to increase the uncertainty in Lboi and 
dilute the connection between observables such as the 
rise time and physical quantities such as the diffusion 
timescale. 

The broad-lined SNe Ic associated with GRBs have 
been extensively studied in the context of aspheri- 
cal models, with different studies coming to confiict- 
ing conclusions for the same object ( e.g., SN 1998bw: 
Hofiich et al. 1999^ IMaeda e t al.' 20061 The models of 
Maeda et al. (2006.) for SN 1998bw predict that lines of 
sight along the major axis of a bipolar outfiow will see 
light curves with faster rise times and higher luminosi- 
ties than average, perhaps indicating a general trend with 
relevance for PSl-lOafx. However, despite the high ex- 
plosion energies (>10^^ erg) and fairly high nickel masses 



(^0.4 M0) assumed in their models, the maximum peak 
luminosity was still an order of magnitude lower than 
for PSl-lOafx. In addition, the spectral features of PSl- 
lOafx are not as blended as those of broad-lined SNe Ic 
(Figure in]). We conclude that an aspherical SN 1998bw- 
like model is unlikely to produce an object that resembles 
PSl-lOafx unless some other ingredient is added to the 
models. 

Another possibility in a CSM interaction scenario is 
that the material with which the SN is interacting is 
distributed in an aspherical fashion, perhaps allowing 
a distant observer to see both the luminous continuum 
emission from the interaction region and the SN ejecta. 
However, such models still require the combination of the 
covering fraction of the CSM, the radiative efficiency fac- 
tor, and the SN kinetic energy to be sufficiently high to 
produce the ~10^^ erg of o ptical eni ission that we mea- 
sure. The SLSN models of lMetzgeri ((2010.) invoke rem- 
nant protostellar disks around massive stars, but do not 
simultaneously fit the high peak luminosity and the fast 
timescale of PSl-lOafx. In addition, such a model has 
similar difficulties as spherical CSM interaction models 
in explaining the strength of the P-Cygni features in the 
observed spectra if there is continuum emission from an 
interaction region that dilutes the strength of the spec- 
tral features from the ejecta. 

6.4. Is PSl-lOafx a Lensed SN? 

We now consider whether gravitational lensing could 
result in an object with the spectrum and light curve 
shape of a normal SN, but artificially boosted in lumi- 
nosity. However, the peak luminosity of PSl-lOafx is a 
factor of >50 higher than normal SNe Ic, requiring an 
extremely high magnification factor despite the lack of 
an obvious lens. While th i s work was being refereed, a 
preprint by iQuimbv et al.l (|2013f ) expanded on this hy- 
pothesis further and proposed that PSl-lOafx was in- 
stead a normal SN la and therefore had to be lensed by 
a magnification factor of ^20. 

There are several constraints against the lensing hy- 
pothesis, starting with the lack of an apparent lens. In 
standard ACDM cosmologies, the optical depth to lens- 
ing with amplification factors >10 is very small, 10^^ 
— 10~^ for sources at the redshift of P Sl-lOafx (e.g., 
iHilbert et aL|[2008l: iTakahashi et al.|[20TTI ). and the mass 
distribution of the lenses for those extreme magnifica- 
tions is peake d at cluster-scale haloes of ~10^^ Mq 
(jHilbert et all HoOS). However, there is no foreground 
galaxy cluster visible in the PSl images, as lQuimbv et al.l 
poll also concluded on the basis of a clustering analysis 
of independent Canada-France Hawaii Telescope images 
of the field. The source present along the line of sight 
to PSl-lOafx exhibits [O II] emission and the SED (in- 
cluding a Balmer break) of a galaxy at the same redshift 
as the SN. The [O II] and rest-frame UV continuum lu- 
minosities produce similar estimates of the SFR. These 
observations are all consistent with the continuum flux 
from the source along the line of sight being dominated 
by, if not entirely due to, emission from the host galaxy 
of PSl-lOafx, and leave little room for a contribution 
from a hypoth etical lens galaxy at l ower redshift. These 
constraints led IQuimbv et al.l ()2013D to propose the exis- 
tence of a dark lens from a previously unknown popula- 
tion, such as a free-floating black hole or a dark matter 
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halo with few baryons. 

This provocative conclusion is based on the premise 
that PSl-lOafx has to be classified as a normal SN la and 
therefore the intrinsic luminosity is known. We demon- 
strated above that the spectra of PSl-lOafx are similar 
to those of both SNe la and Ic near maximum light (Sec- 
tion 4) . We also note that the presence of an absorption 
feature from Si II A6355 does not uniquely classify an ob- 
ject as a SN la because that line is also present in some 
SNe Ic. As an experiment, we took near-maximum light 
spectra of an unambiguo us nearby SNe Ic tha t exhibits 
Si II features, SN 2004aw (jTaubenberger et al . 2006). re- 
sampled the data, and added noise to match the S/N of 
our PSl-lOafx FIRE data that show the possible absorp- 
tion feature near 6100 A. These noisy spectra were then 
run th rough the SuperNov a IDentification (SNID) code 
of B londin fc Tonrvl (l200l , using the updated templates 
of [Silverman et al.l (|2012[ ). with care taken to exclude 
templates from the same object. We found that typi- 
cally all of the top 20 best matching templates are SNe 
la. This is in part due to the biased sample of spectral 
templates used by SN classification codes, where SNe la 
a re largely overreprese nted. 

iQuimbv et al.l ()2013[ ) also claim that the light curve of 
PSl-lOafx is a good match for SNe la. However, we note 
that in order to get a good reduced-x^ for their fit, they 
had to artificially increase the error bars in the NIR and 
exclude the J-band data. Furthermore, all of the PSl- 
lOafx observations prior to a phase of —5 d fall below 
the SN la prediction, including the zpi non-detection at 
—7.4 d that is ~0.7 mag below their model. This is a 
consequence of the fast rise time of PSl-lOafx, which we 
measure above to be 11.8 ± 0.7 d if we assume that the 
rise scales as Lboi oc and 8.7 ± 1.9 d if we relax that 
constraint. Although spectroscopically normal SNe la 
show a small dispersion in their rise times, objects with 
normal light curve shapes have rise times of 16—18 d in 
B (e.K.,|Hav den et al. 2010; GancshaUnEam et al.ll2011l: 
IGonzgilez-Gaitan et al.ll2012i ). iJha et al.i (\200W demon- 
strated that the rise in U (closer in rest wavelength to our 
zpi observations) was 2.3 d shorter, but still significantly 
longer than observed for PSl-lOafx. At phases of —12 to 
— 10 d , normal SNe la emit fluxes about 20—40% of the 
peak (jGaneshalingam et al.|[20Tl T). but that is about Aa 
above our non-detections (Figure [T^ . 

In summary, although SNe la have reasonably similar 
spectra and light curve shapes to PSl-lOafx, the details 
of the rise time are not a good fit for SNe la and the spec- 
tra of the unambiguous SN Ic 2004aw are also similar to 
SNe la when degraded to our low S/N ratio. Therefore, 
an identification of PSl-lOafx with normal SNe la is not 
required by the data, and there is no need to invoke a 
new population of unseen dark gravitational lenses. In- 
stead, we accept the lack of an apparent lens at face value 
and conclude that PSl-lOafx was actually an unusually 
luminous hydrogen-deficient SN. 

7. CONCLUSIONS 

We present multiwavelength observations of PSl-lOafx 
at redshift z = 1.388, perhaps the most luminous 
SN yet discovered. The combination of observables 
presents strict constraints on any theoretical interpre- 
tation. These are: 
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• A peak bolometric luminosity of 4.1x10'*'' erg s~^ 

• An observed rise time of '^12 d 

• A fast light curve decay, with Ami^— 0.95 mag 

• At least 7x10^° erg of optical radiation emitted, 
with a total likely closer to 10^^ erg 

• Red color near maximum light (Tbb = 6800 K), 
although Tcs may be closer to 12, 000 K 

• Constant UV color before maximum light 

• Photospheric velocities near maximum light of 
-11,000 km S-* 

• Spectra that most closely resemble normal SNe Ic, 
with deep Ca II P-Cygni absorption 

• Photometric and spectroscopic evidence for some 
line blanketing in the NUV from iron-peak ele- 
ments, although not as much as for SNe la 

• A massive (—2x10*° M©) host galaxy that is un- 
likely to have an extremely low metallicity 

We surveyed existing models for SLSNe and found 
that none were acceptable. In particular, the large in- 
ferred i?BB near maximum light is very hard to reconcile 
with the fast observed rise time and the measured pho- 
tospheric velocities because the SN ejecta need too much 
time to reach such large radii. The magnetar models that 
match the peak luminosity and rise time do so by pro- 
ducing a much higher temperature at a smaller radius. 
These fiaws are generic to models powered by internal 
energy sources, although if the onset of the internal en- 
ergy source can be sufficiently delayed after the initial 
explosion, then maybe the conflict between the veloci- 
ties and radii can be avoided. Shock breakout scenarios 
invoking dense CSM provide a promising solution by al- 
lowing the SN to explode well before the optical emis- 
sion is detectable. If the CSM is sufficiently dense out 
to —5x10*^ cm, then the light curve will only rise after 
the ejecta have had time to reach that radius. However, 
these models leave the normal SN Ic-like spectrum un- 
explained and are in apparent conflict with the observed 
color evolution before maximum light. More detailed ra- 
diative transfer calculations of the models are necessary 
to know whether these fiaws can be avoided. In addition, 
it is not clear how to produce the special CSM structure 
that has to be assumed (a truncated wind or a shell). 

PSl-lOafx was initially identified as an object of in- 
terest due to its unusual colors (rpi — zpi w 2 mag at 
peak) and would not have been found without the PSl 
observation strategy, which includes regular ipi and zpi 
observations as part of the search. At redshifts below 
— 1, an object like PSl-lOafx would have a strong de- 
tection in rpi and would not stand out without knowing 
the redshift. Unlike many of the SN 2005ap-like SNe, the 
time-dilated light curve of a PSl-lOafx-likc SN at lower 
redshift would not be unusually long and would not be 
associated with the attention-grabbing lack of a visible 
host galaxy. However, PSl-lOafx was more than a mag- 
nitude brighter than its host galaxy, which potentially 
provides a selection criterion to increase the odds of find- 
ing similar objects in the future. In the sample of SLSNe 
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found by PSl, SN 2005ap-like objects (iChomiuk et alj 
l201lHBerger et al.ll2012HLunnan et al.ll2013D outnumber 
PSl-lOafx by at least an order of magnitude. Our poor 
understanding of spectroscopic incompleteness and the 
relevant selection effects precludes a more precise esti- 
mate of rates at this time. 

An additional complication to be considered is the 
role of metallicity. The massive host of PSl-lOafx 
is in contrast wit h the low-niass (iNeill et al.l [20111 
and low-metal l icitv (lYoung et al.l[2010l : lChen et al. I I2OI 
iLunnan et al.l 120131 ) hosts of known hydrogen-poor 
SLSNe. The natal composition of the progenitor could 
have an indirect effect on the explosion through stellar 
evolutionary processes or it could directly affect the ap- 
pearance through the opacity in the outer ejecta (could 
SN 2005ap-like events with higher metal abundances ex- 
hibit line-blanketed spectra like PSl-lOafx?). 

Most of the published SLSNe that lack hydrogen in 
their spectra are siniilar to SN 2005ap and SCP06F6 
(iQuimbv et al.ll2007t [Barbarv et al.112 009': 'Quimbv et al " 



20111: iPastoreho et al.l 120 lOt ICho miuk ct al. 2011; 



Leloudas et al .l [20121: ILunnan et al.' 2013^. SN 2007bi 
^al-Yam ct atl 120091: [Young ct al. 2010) and a couple 
of as-yet unpublished objects (|Gal-'Yainll2012[ ) were the 
only known exceptions. They exhibited rather different 
spectra and had ver y slow light cu r ve de cays, which 
were interpreted by (Cal-Yam etall (|2009D to be the 
result of decay of a large amount of radioactive ^^Ni 
produced in a pair-instability explosion. This result has 
been ch allenged on theoretical grounds bv lDessart et al.l 
(|2012b[) , who show that the observations do not match 
theoretical expectations for pair-instability SNe. With 
the discovery of PSl-lOafx, we now have further evidence 
that the most luminous SNe are a heterogeneous lot, 
even when just considering the objects lacking hydrogen. 
A key question for future investigations is whether these 
diverse outcomes of stellar evolution can be produced 
by variations on a single underlying physical model 
or whether multiple pathways exist to produce SNe of 
these extraordinary luminosities. PSl-lOafx differs from 
existing SLSNe-I in almost every observable except for 
the peak luminosity. It is hard to understand how all 
of these differences could be produced with only small 
modifications to a single model, and may indicate that 
PSl-lOafx is the first example of a different channel for 
producing SLSNe. 
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